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I.  INTRODUCTION 


This  is  the  third  g^arte^ly  report  issued  under  Contract 
N00017-71-C-1417.  Work  under  the  contract^ is  aimed  at  determining 
the  feasibility  of  implementing  adaptive  echo  cancellation 
techniques  to  underwater  weapons.^ 

In  work  reported  in  progress  reports  1 and  2,  effort 
was  aimed  at  determining  the  dynamic  characteristics  of  the  adaptive 
system  represented  in  Fig.  1.  In  its  major  attributes  thiB  system 
is  identical  to  one  utilized  by  M.  H.  Sondhi1  to  cancel  "echoes" 
generated  at  discontinuities  in  telephone  communications  lines. 

A more  complex  form  of  the  system  haB  been  adopted  for  study 
within  the  last  quarter,  but  the  system  described  in  Fig.  1 is 
a convenient  vehicle  for  introductory  remarks  concerning  the 
system's  operation. 

The  function  F(t)  represents  a received  echo  signal 
that  is  representible  as  a manifold  of  time  shifted  components. 

The  components,  in  turn,  are  regarded  as  time  shifted  replicas 
of  the  reference  function  R(t) . 

An  examination  of  Fig.  1 will  show  that  the  reference 
signal  R(t)  is  admitted  into  a tapped  delay  line  with  uniform 
delay  increment  t between  consecutive  taps.  The  output  signal  at 
each  tap  is  therefore  a time  shifted  replica  of  the  transmit 
reference  function  and  this  signal  is  connected  as  one  input  to 
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Sondhi,  M.H.,  An  Adaptive  Echo  Canceller,  B.S.T.J.  46, 
March  1967,  pp.  497-511 
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a multiplier.  The  second  input  to  the  multiplier  is  an  error 
signal  e(t),  that  is  common  to  all  the  multipliers  in  the  system. 
The  output  of  any  one  multiplier  forms  the  input  to  an  integrator, 
and  an  examination  of  Fig.  1 will  show  that  the  multiplier- 
integrator  components  implement  a cross- correlation  between  the 
delayed  reference  signals  R(t-nT)  and  the  common  error  signal 
e(t).  The  output  of  the  integrator  is  admitted  aB  a control 
voltage  to  a second  multiplier  unit  gn>  The  delayed  reference 
signal  is  connected  as  a second  input  to  this  multiplier  and  the 
output  of  the  unit  can  be  expressed  as  g^R(t-n*r).  The  gain 
parameters  gn  can  take  on  positive  or  negative  values  and  can 
also  take  on  the  value  zero.  As  is  indicated  in  Fig.  1,  the  value 
of  gn  is  controlled  by  its  associated  cross-correlator.  The 
output  signals  of  the  correlator  units  are  summed  together  and 
this  composite  signal  is  subtracted  from  the  received  signal, 

F(t),  to  form  the  error  signal  e(t).  Thus,  under  ideal  cir- 
cumstances, the  correlators  provide  a measure  of  the  time 
registry  between  their  reference  voltage  outputs  and  the  com- 
ponents of  the  echo  signal  and  control  the  gain  parameters  to 
null  the  error  signal.  The  study  of  the  system's  response 
under  less  than  ideal  conditions  is  a main  goal  of  the  work  in 
progress. 

The  intent  is  to  use  the  processor  to  implement  what 
may  be  termed  as  a steady  state  mode  of  operation.  This  mode 
of  operation  requires  that  each  of  the  signals  taken  from  the 
tapped  delay  line  be  in  an  active  state  over  the  whole  time 
interval  of  a given  test.  This  requirement  imposes  constraints 
on  the  time  relationship  between  F(t)  and  the  reference  signals, 
and  also  stipulates  a maximum  time  duration  for  a given  test. 

The  onset  of  F(t)  will  define  the  zero  of  time  as  is 
indicated  in  Fig.  2.  The  requirement  that  all  of  the  delayed 
reference  signals  be  in  an  active  state  during  the  time  span 
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of  a test  implies  that  the  reference  functions  be  in  an  active 
state  before  the  onset  of  F(t) . Moreover,  examination  of  Fig.  2 
will  show  that  the  time  interval  over  which  steady  state  tests 
can  be  conducted  is  specified  by 


0 < t * T , - T . 

a s 


By  inspection  of  Fig.  1,  the  error  signal  can  be 
described  by 


e(t)  » F<t)  - 2 9k(t)Rk(t)  , 
k®0 


where 


R^t)  » R(t~kT) 


In  addition,  an  examination  of  Fig.  1 will  show  that  for  each 
of  the  integrator-multiplier  units  one  can  write  an  equation 
of  the  form 


i rt 

9n(t)  “ gn(0)  + T J Rn(xMx)dx 


The  Bystem's  feedback  is  made  more  apparent  by  taking 
the  time  derivative  of  Eq.  (4)  and  writing  out  the  full  ex- 
pression for  the  error  function  e(t).  This  results  in  N 
equations  of  the  form 


imLSUiciH 


M"*4. 

o R_(t) 

9n(t)  - -^CFtt)  - 2u  SkV 


In  the  last  progress  report,  numerical  techniques  were 
described  that  permitted  a direct  integration  of  the  coupled 
equations  represented  by  Eq.  (5).  The  computer  program,  as 
described  in  the  second  progress  report  entailed  the  following 
operations: 

1.  A first  phase  of  the  program  was  concerned 
with  obtaining  accurate  values  of  the  gain 
parameters  at  four  aonsecutive  points  in 
time  through  a process  of  iteration, 

2.  A second  phase  of  the  program  was  con- 
cerned with  continuing  the  solution  through 
a process  of  extrapolation. 

Early  exercises  of  the  program  clearly  indicated  that 
the  extrapolation  technique  did  not  provide  stable  solutions  to 
the  coupled  equations  represented  by  Eq.  (5) . To  remedy  this 
defect  an  additional  iterative  operation  was  required  after 
step  2,  as  described  above.  With  this  modification  the  program 
operated  satisfactorily  and  was  used  to  demonstrate  to  this 
program's  sponsor  the  computed  response  of  the  system  for  a 
number  of  cases.  These  cases  are  touched  on  briefly  in  the 
following  section  of  the  report.  However,  an  additional 
deficiency  became  apparent,  this  time  relating  to  the  configura- 
tion of  the  system  described  in  Fig,  1.  The  system,  as  described, 
is  adequate  only  for  the  investigation  of  those  cases  in  which 
the  signal  components  are  in  exact  time  registry  with  a subset 
of  the  reference  signals  taken  from  the  tapped  delay  line. 
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Accordingly,  agreement  was  reached  with  the  program  sponsor  to 
generalize  the  configuration  of  the  system  through  the  adoption 
of  a Quadrature  Loop  system.  Description  of  this  system  is 
deferred  to  a later  section  of  this  report.  However,  investiga- 
tion relating  to  the  new  system  finally  entailed  the  development 
of  a new  computer  program.  This  program  is  again  based  on 
iterative  techniques,  but  differs  substantially  in  detail  from 
the  original  program  developed  for  use  with  the  system  represented 
in  Pig.  1.  Moreover,  through  slight  modification,  the  new  program 
can  be  used  for  investigation  of  either  the  old  or  the  new  system 
configurations.  Accordingly,  further  discussion  of  the  program 
described  in  the  second  progress  report  will  be  omitted.  Finally, 
results  computed  with  the  new  program  are  presented  and  discussed 
in  the  following  pages  of  the  report. 


II.  TECHNICAL  PROGRAM 
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RE SULTS  COMPUTED  ON  THE  BASIS  OF  SONDHI's  SYSTEM 


Before  addressing  the  Quadrature  Loop  system,  a number 
of  results  computed  for  the  system  configuration  shown  in  Fig.  1 
will  be  presented.  The  results  in  question  are  solutions  of  the 
coupled  gain  equations  represented  by  Eq.  (5)  and  were  obtained 
through  use  of  a modified  form  of  the  computer  program  described 
in  the  second  progress  report. 

Figure  3 is  a plot  of  the  gain  versus  time  variation 
for  a single  loop  system  in  which  the  received  signal  is  taken 
as  being  identical  to  the  reference  signal  and  with  both  signals 
occurring  in  perfect  time  registry.  The  case  in  point  is  useful 
in  demonstrating  that  each  of  the  system  loops  can  be  regarded 
as  having  the  frequency  response  characteristics  of  a single  pole, 
low-pass  filter. 

For  a single  loop,  Eqs.  (5)  take  the  form 


g0  “ - g0V  * 


(6) 


Through  a rearrangement  of  terms  and  the  introduction  of  an 
integrating  factor,  Eq.  ( 6 ) can  be  written  as 


m 


Fig.  3. 


it  15  to  25  30” 

Time  (msec) 


Gain  Versus  Time  for  Single  Loop,  Sondhi  Type 
System,  Responding  to  Single  Component  Signal 
of  Zero  Doppler  with  Perfect  Time  Registry. 
Received  and  Reference  Signals  are  Linear  FM. 
Slides  of  0.064  Sec.  Duration,  1 kHz  Bandwidth, 
Start  Frequency  of  0.5  kHz.  Integrator  Time 
Constant  is  4 msec. 
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Equation  ( 7 ) can  then  be  integrated  to  yield 


J Ro(^)dTi 


"t  I Ro(,n)dTi 


g0(t)  - g0(o) 


+ £ J R0(x)F(x)  e x 


Because  RQ(x)  is  an  PM  slide,  it  can  be  represented  in  the  form 


Rq(x)  - cos  ©0(t)  , 

where 

9q (t ) » U)Qt  + ^t2 


In  Eq.  (10),  ujq  represents  the  beginning  radian  frequency  of  the 
FM  slide  and  p is  a measure  of  the  signal's  bandwidth.  Because 
of  the  trignometric  form  of  the  right  hand  side  of  Eq.  ( 9 ) , the 
square  of  the  function  can  be  expressed  as 


Rq  (ll>  a | | COS  290(n)  . 


The  second  term  on  the  right  hand  of  Eq.  (l 1 ) represents  the 
second  harmonic  of  the  PM  slide.  Because  of  the  oscillatory 
nature  of  this  term  it  will  be  assumed  that  the  time  integral 

p 

of  Rn(r|)  can  be  approximated  by  the  integral  of  the  constant  term, 
1 ^ 

•jg,  and  that  the  integral  of  the  second  term  on  the  right  in 
Eq.  (11)  can  be  ignored.  With  this  approximation  Eq.  (8)  can 
be  written  in  the  form 


— ■■■■■■■■■; " 1 ''  ''  ■■  ■'  ..  ^.ll: jiitVtiy,: is; 


g0(t)  90(0> 


“ 2T  1 **  2T1  ^**^) 

e 2T  + f J e 2T  RQ(x)F(x)dx  . 

0 

(12) 

Examination  of  the  second  member  on  the  right  in  Eq.  (12)  will 
show  that  the  term  has  the  form  of  a convolution  integral.  The 

same  expression  describes  the  output  of  a filter  that  has  an 

« X 

2T 

impulse  response,  e , and  an  input  signal  equal  to  the  product 
Rq(x)F(x)  along  with  an  impulse  function  to  generate  the  term 
g(0).  A filter  of  this  type  can  be  realized  with  a series  con- 
nected RC  network  with  the  input  signal  applied  to  the  resistor 
and  the  output  signal  taken  across  the  capacitor.  The  realiza- 
tion requires  only  that  the  RC  product  be  equal  to  2T.  Thus, 

if  the  filter's  low  frequency  cut-off  is  designated  by  f , then 

c 


f 


c 


2rrRC 


4tiT 


(13) 


If  F(x)  is  taken  equal  to  R(x),  and  if  the  product  is  again  ap- 
proximated by  the  factor  -ip  then  substitution  in  Eq.  (12)  yields 

X-  t 

90(t)  « gQ(0)  e 2T  + [1  - e 2T]  . (14) 


As  a test  of  the  approximations  leading  to  Eq.  (14),  the  computer 
program  was  used  to  compute  gQ(t)  for  a case  in  which  gQ(0)  was 
zero,  F (t)  was  equal  to  RQ(t),  with  Rg(t)  described  by  Eqs.  (9) 
and  (10).  In  addition,  the  beginning  frequency  of  the  FM  slide 
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was  taken  at  500  Hz,  and  the  signal's  bandwidth  was  taken  as 
1000  Hz  over  a 64  msec  duration.  The  integrator  time  constant 
T was  set  equal  to  4 msec,  and  the  computed  results  are  presented 
in  graphical  form  in  Pig.  3 . The  results  obtained  with  the 
computer  are  plotted  as  a solid  curve  while  points  evaluated 
from  Eq.  (14)  are  represented  by  the  symbol  x along  the  curve. 

The  scale  to  which  the  graph  is  drawn  does  not  permit  a visual 
distinction  between  the  values  arrived  at  by  the  two  methods 
described,  and  a more  careful  examination  of  numerical  data  in- 
dicates that  the  deviation  between  the  results  obtained  with 
the  computer  versus  those  computed  from  Eq.  (14),  is  on  the  order 
of  0.05  percent.  Thus,  the  approximations  and  conclusions  en- 
tailed in  the  development  of  Eq.  (14)  are  considered  valid. 

1 

The  computer  program  developed  for  the  investigation 
of  Sondhi's  system  was  used  to  explore  a number  of  cases,  and 
certain  of  those  cases  were  run  with  the  project  sponsor  in 
attendance.  Because  current  interest  relates  to  a modified  form 
of  the  system,  computed  results  for  only  one  additional  case 
involving  Sondhi's  system  will  be  presented  in  this  report.  The 
case  involves  a four-loop  system  responding  to  a three-component 
signal.  The  components  were  considered  to  have  levels  that  were 
respectively  8,  4,  and  2 times  the  levels  of  the  system's  reference 
functions.  Moreover,  the  components  were  considered  to  be  un- 
altered by  Doppler  and  occurred  in  perfect  time  registry  with 
the  first  three  corresponding  reference  signals  of  the  system. 

The  computed  results  are  presented  in  graphical  form  in  Fig.  4. 

The  main  point  in  presenting  these  data  is  to  show  that  mutual 
interference  occurs  between  the  loops  of  Sondhi's  system.  This 
action  is  evidenced  by  a high  frequency  ripple  superimposed  on 
a slowly  varying  component  associated  with  each  of  the  gain  para- 
meters. The  high  frequency  component  is  evident  in  the  drawing's 


of  Fig.  4 , and  is  attributed  to  a given  loop's  response  to  a 
signal  component  not  in  time  registry  with  the  reference  function 
associated  with  the  loop.  However,  Fig.  4 shows  that  for  the 
case  considered,  the  system  still  proceeds  to  a null  solution 
although  not  in  a monotonic  manner.  A similar  phenomenon  is 
observed  with  a modified  system  described  later  in  this  report, 
and  to  avoid  repetition,  further  discussion  of  this  phenomenon 
will  be  deferred  pending  description  of  the  newer  system. 
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1 1 /kmseo/l  Reference  Onset  Distribution 


8j 4J! 2 '/[ Signal  Component  Onset  Distribution 


Time  (msec) 

Fig.  4.  Gain  Parameters  Versus  Time  for  Four  Loop,  Sondhi  Type  Syst 
Responding  to  Three  Component  Signal  with  8,  4,  2 Amplitude 
Distribution,  Data  Shown  Are  for  Zero  Doppler  and  Perfect 
Time  Registry.  System  and  Signal  Parameters  Are  as  Cited  i 
Fig.  3. 


15 


r 


M 


| 

fct 


I 


B.  THE  QUADRATURE  LOOP  SYSTEM 

Owing  to  the  fact  that  the  echo  canceller's  error 
signal  nulls  through  a subtractive  process,  a precise  align- 
ment  in  time  must  be  maintained  between  a subset  of  the  ref- 
erence signals  and  the  various  components  of  the  echo  in  order 
to  achieve  a nulling  action  in  Sondhi's  system.  Thus,  in 
practice,  misalignments  in  time  registry  on  the  order  of  one 
fourth  period  at  the  center  frequency  of  the  FM  signals  would 
be  enough  to  prevent  the  nulling  action,  because  this  situation 
would  be  roughly  equivalent  to  subtracting  two  sinusoids  charac- 
terized by  a relative  phase  shift  - a process  which  also  cannot 
produce  a null  that  is  uniform  with  time.  Accordingly,  through- 
out the  last  quarter,  studies  have  been  baaed  on  a modified 
form  of  Sondhi's  system  that  incorporates  new  degrees  of  freedom 
to  enable  an  adaptive  form  of  phase  compensation.  The  modified 
echo  canceller  is  referred  to  as  a Quadrature  Loop  system  through- 
out this  report  and  its  configuration  is  described  by  Fig.  5. 

As  can  be  inferred  from  a comparison  of  Figs.  5 and  1 
the  quadrature  loop  system  differs  from  Sondhi'B  in  that  two 
reference  functions  are  supplied  to  the  system  and  each  loop  is 
provided  with  two  separate  paths.  The  reference  function  R(t) 
is  expressed  as  before  by 


R(t)  » cos(u>0t  + ti|-) 


(15) 


H \ 
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while  the  function  X(t)  is  the  quadrature  component 


X(t)  = sin(wQt  + . 


(16) 


Tapped  Delay  Lines 


Adaptive 


The  two  functions  are  represented  in  Fig.  5 as  being  transmitted 
along  separate  tapped  delay  lines.  The  signals  taken  after  n 
delay  elements  can  therefore  be  represented  as 


Rn(t)  - R(t-nT)  , and  (17) 

Xn(t)  ■ X(t-nt)  . (18) 


At  each  such  location  two  adaptive  gain  control  units 
similar  to  those  employed  in  Sondhi's  system  scale  the  signal 
levels  transmitted  through  the  loops.  Thus,  the  contribution  to 
a central  summation  point  made  by  the  quadrature  loop  energized 
after  n delay  elements  can  be  expressed  as 


Vk>  ■ + WC> 


(19) 


The  mechanism  through  which  adjustment  of  the  gains 
an#  Pn  can  result  in  an  adjustment  of  phase  angle  can  be  most 
readily  demonstrated  by  first  expressing  the  gain  parameters  in 
terms  of  auxiliary  quantities  Cn,  $n,  through  the  equations 

an(t)  - Cn(t)  cos  *n(t)  (20) 

pn(t)  - Cn(t)  sin  »n(t)  . (21) 


Then  expansion  of  Eq.  (19)  yields 
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sn(t)  “ Cn(t)[cos  *n(t)  cos[w0(t-nT)  + ^ (t-m-)2] 


+sin  Sn(t)  sin[iw0(t~nT)  + ^ (t-nt)2]} 


(22) 


Equation  (22)  can  be  recognized  as  being  trignometrically 
identical  to 


Sn(t)  ■ Cn(t)  cos[oJ0(t-nT)  (t-m-)2  - *n(t)]  . 

(23) 

Equation  (23)  of;  course,  provides  no  information  concerning  the 
system's  adaptive  response.  Nevertheless  , considerable  insight  i 

can  be  gained  with  respect  to  the  computer  generated  data  presented 
in  a later  section  of  this  report  by  first  considering  a few 
elementary  but  highly  hypothetical  cases. 

As  a first  case,  consideration  will  be  directed  toward 
a single  quadrature  loop  system  responding  perfectly  to  a single 
component  echo  but  with  the  echo  component  misaligned  in  time 
with  the  reference  signals  supplied  by  the  tapped  delay  line. 


i- 


For  convenience,  take  the  quadrature  loop  signal  as 


SQ(t)  » CQ(t)  cos[u)0t  + ^t2  - *Q(t)]  , (24) 


and  take  the  echo  signal  in  the  form 
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F (t)  * cos[u>0(t-e)  + ^(t-e)2} 


cos[(u)nt  + ¥t2)  - (net  + - ^e2)]  . 

0 2 0 2 (25) 


The  second  expression  on  the  right  hand  side  of  Eq.  (25)  is  ar- 
rived at  by  a simple  rearrangement  of  terms.  However,  comparison 
of  this  expression  with  Eq.  (24)  will  show  that  the  two  can  be 
made  identical  by  setting 


*0(t) 


net  + u!0e 


Cn(t)  «*  1 


Thus,  by  Eqs.  (20)  and  (21),  the  gain  parameters  for  the  loop 
would  be 


an  « cos(|iet  + ujQe  - ^e2) 


5n  ■ sin(^et  + u)Qe  - ^e2) 


Thus,  in  a perfect  system  of  the  type  considered,  the  gain  para- 
meters would  describe  a sinusoidal  variation  with  time  with  a 
radian  frequency  equal  to  ne.  One  can  therefore  see  that  the 
frequency  of  the  variation  increases  directly  with  the  time  mis- 
alignment e . 
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As  a second  case,  consider  a single  component  echo 
with  Doppler  that  is  in  exact  time  registry  with  a single  quadrature 
loop  hut  with  reference  functions  that  do  not  reflect  Doppler. 

For  this  case  take 


SQ(t)  * CQ(t)  cos [ uoQt  + ^t2  - »0(t>]  , 


and  assume  the  received  signal  given  by 


F(t)  >=.  cos[au)0t  + na  ' 3 


In  Eq.  (31)  the  ef feats  of  Doppler  are  reflected  by  applying  a 
constant  scale  factor  a to  the  time  parameter.  By  equating  the 
arguments  of  the  right  hand  functions  of  Eqs.  (30)  and  (31),  one 
can  readily  show  that  the  two  equations  are  identical  if 


*0(t)  «*  (u0(l-a)t  + ^(l-a2)t2 


Cft(t)  « 1 


Equation  (32)  shows  that  the  gain  parameters  for  this  case  would 
describe  a linear  FM  variation  with  time. 

As  noted  previously,  the  two  cases  juBt  considered  do 
not  actually  describe  the  dynamic  response  of  the  system.  A com- 
puter program  has  been  utilized  to  demonstrate  system  response 


for  a number  of  cases.  Referring  to  Fig.  5 , the  gain  para- 
meter  in  any  loop  is  the  integral  of  the  product  of  the  corres- 
ponding reference  signal  and  the  error  signal.  Thus,  the  gain 
parameters  as  a function  of  time  can  be  inferred  by  solving  the 
coupled  linear  differential  equations  specified  by 

-2cVk  + W1  ' <34> 

k 

K " • <33> 


I 

• t 


N • 


The  computer  program  in  question  implements  a direct  numerical 
integration  of  the  coupled  equations  represented  by  Eqs.  (34) 
and  (35) . 


Kt ...  f 


IUinIitejclH 


C.  BRIEF  DESCRIPTION  OF  THE  COMPUTER  PROGRAM  USED  FOR 


INVESTIGATION  OF  THE  QUADRATURE  LOOP  SYSTE 


The  computer  program  developed  for  use  in  an  investiga- 
tion of  the  Quadrature  Loop  system  resembles,  in  many  respects, 
the  program  developed  for  a study  of  Sondhi's  system.  This  being 
the  case,  a detailed  description  of  the  new  program  will  be 
deferred  to  the  final  report.  The  salient  features  of  the  new 
program  are  outlined  below. 

The  operation  of  .the  program  is  based  on  numerical 
methods  of  solution  derived  in  the  second  progress  report.  Only 
the  results  of  those  derivations  will  be  presented  here,  and  are 
expressed  by  the  following  equations: 


am(Y+h)  « am(Y)  + {5Rm(v>e(Y>  + BRm( y+h)  e (Y+h) 


-Rm(Y+2h)e(Y+2h)} 


Bm(Y+h)  - Pm(Y)  + iff  (5Xm(Y)e(Y)  + 8Xm ( Y+h) e ( Y+h) 


-Xm(Y+2h)e(Y+2h)} 


am(Y+2h) 


cttn(Y)  + [Rm(Y)e(Y)  + 4Rm(Y+h)e(Y+h) 


+Rm(Y+2h) e ( Y+2h) } 
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Sm(V+2h)  ■=  Um(y)  + j~  [Xm(y>  eCv)  + 45£m(Y+h)c(y+h) 


+Xm(Y+2h) c(v+2h) ) . (39) 


In  Eqs.  (36)  through  (39),  a^,  0m  are  the  gain  parameters,  Rm# 
Xm  are  the  reference  functions,  and  T is  the  integrator  time 
constant.  Moreover,  ^ is  a particular  value  of  time  and  h is 
a small  time  increment.  The  error  function  is  given  by 


® (Y)  - F(y) 


-E 


Cak(y)Rk(Y)  + Pk(Y)Xk(y)] 


(40) 


where  F(y)  is  the  received  signal  at  t-y. 

Input  to  the  program  includes  initial  values  for  the 
gain  parameters  at  t*=0.  The  first  phase  of  the  program  is  aimed 
at  finding  values  of  the  gain  parameters  at  t»=h  and  t=2h.  This 
is  carried  out  through  a process  of  iteration  involving  Eqs.  (36) 
through  (39)  with  y set  equal  to  zero.  On  the  right  hand  sides 
of  these  equations,  the  iterations  relate  to  the  terms  e(h)  and 
c(2h).  The  process  commences  by  providing  initial  estimates  for 
ak(h),  ak(2h),  0k(h),  and  pk(2h).  Then  through  use  of  Eq.  (40) 
initial  estimates  of  e (h)  and  e (2h)  are  computed.  Iteration  then 
commences  by  computing  a new  value  for  a0(h)  from  Eq.  (36)  after 
which  the  initial  estimate  of  e (h)  is  corrected  to  reflect  the 
new  value  of  aQ(h).  The  same  process  continues  for  0Q(h),  a0(2h), 
P0(2h)  and  commences  for  a^(h).  One  iteration  is  complete  after 
the  process  has  propagated  through  0N>>i^(2h),  where  N is  the  number 
of  Quadrature  Loops,  A test  is  then  applied  to  determine  the 
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error  between  successive  iterations.  Thus,  if  n in  onth) 

II  Jv 

designates  the  n iteration,  the  error  is  specified  by 


ErII  ln+lak<mh)  - nak(mh)l  + ln+l®k(n,h)  - nMmh)l  • 
H,=1  (41) 

During  each  iteration  a new  value  for  is  computed  and  the 

iterative  process  continues  until  E1  is  less  than  a maximum  value, 

-10  ^ 

typically  taken  as  10  . In  most  cases  the  error  falls  below 

this  maximum  during  five  iterations  carried  out  for  a four  loop 
system. 

After  refined  values  become  available  for  the  gain 
parameters  at  t=0,  t=h,  and  t=2h,  the  program  implements  an 
iterative  process  based  on  Eqs.  (38)  and  (39)  alone.  On  the 
right  hand  side  of  these  equations,  the  iterative  process  relates 
solely  to  the  term  e (y+2h) , The  process  is  again  initiated  by 
providing  an  estimate  for  the  gain  parameters  at  the  time  y+2h. 

The  convergence  rate  of  the  iterative  process  is  greatly  improved 
by  solving  Eqs.  (38)  and  (39)  explicitly  for  om(v+2h),  flm(Y+2h) 
on  the  assumption  that  parameters  of  index  other  than  m have  the 
correct  value.  For  example,  in  Eq.  (38),  this  consists  of  solving 

£or  n+l“m<Y+2h)  in 


+ 3T  Rm*Y+2h^n+lait^Y+21’* 


“ am^Y)  + 3T  CRm < V>  ® ( Y)  + Rm(Y+h)e (Y+h) 


+ Rra(Y+2h)em(Y+2h)] 
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where 

em(v+2h)  «c(Y+2h)  + nam(Y+2h)I^(Y+2h)  . (43) 


The  process  represented  by  Eqs.  (42)  and  (43)  acts  to  remove  the 
dependence  on  a^Y*2*1)  from  the  right  hand  side  of  Eq.  (32)  . 

The  reference  functions  and  the  received  signal  functions 
are  generated  within  the  program  through  use  of  signal  generator 
subroutines.  In  cases  involving  Doppler,  the  time  parameter 
within  these  subroutines  is  appropriately  scaled  to  reflect 
assumed  radial  speed  components. 

In  the  cases  considered  in  the  following  section  of 
the  report  computation  was  carried  out  through  use  of  a Data 
General  minicomputer  incorporating  a software  "multiply-divide" . 
subroutine  and  reading  out  via  a standard  teletype  unit.  The 
starting  frequency  of  the  FM  signals  was  taken  to  be  750  Hz 
and  the  frequency  was  assummed  to  slide  through  500  Hz  in  a 
128  msec  time  interval.  For  these  computations  the  increment  h 
was  taken  to  be  50  microseconds,  and  experience  indicates  that 
extending  the  range  of  the  integrations  over  one  millisecond 
requires  approximately  40  seconds  of  computation  followed  by 
about  20  seconds  of  readout.  Therefore  a typical  run  requires 
in  excess  of  1*3  hours  for  execution. 


,D.  COMPUTED  RESPONSE  OF  A FOUR  CHANNEL,  QUADRATURE  LOOP 

SYSTEM  FOR  SEVERAL  TYPES  OF  INPUT  SIGNALS 

1.  Parameters  of  The  System 

The  computed  results  presented  in  the  following 
pages  of  this  report  were  based  on  a system  incorporating  four 
quadrature  loops  and  three  pairs  of  delay  elements.  Throughout 
the  runs  the  time  delays  brought  about  by  the  delay  elements  was 
assumed  to  bo  10  msec.  The  integrator  time  constant  T was  taken 
to  be  4 msec. 

The  FM  signals  used  in  the  computations  were 
assumed  to  commence  at  an  instantaneous  frequency  of  750  Hz 
and  to  increase  linearly  with  time  up  to  1250  Hz  in  a time  in- 
terval of  128  msec.  Thus,  the  signals  were  assumed  to  vary  over 
a 500  Hz  band  in  128  msec.  The  integration  interval*  h,  was 
set  equal  to  50  nsec  throughout  the  computation.  In  all  cases 
the  peak  level  of  the  reference  signals  were  assumed  to  establish 
a unit  level,  and  the  levels  cited  for  the  components  of  received 
signals  are  their  peak  levels  specified  in  terms  of  unit  levels. 

The  data  are  presented  in  the  form  of  manually 
plotted  graphs  of  the  gain  parameters  versus  time.  The  graphs 
were  drawn  from  a data  printout  executed  at  1 msec  intervals. 
Conversion  of  the  printed  data  into  graphical  form  proved  a 
very  laborious  nrocedure  and  this  fact  limited  the  number  of 
cases  that  could  be  included  in  this  report. 

Finally,  a comment  concerning  the  components  of  the 
received  signal  is  in  order.  Each  component  is  considered  made 
up  only  of  a time  shifted  cosine  dependence  such  as  characterizes 
the  reference  functions  :«k(t).  Consequently,  in  cases  involving 
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exact  time  registry  between  undoppled  signal  components  and 
aertain  of  the  reference  signals,  an  exact  null  condition 

requires  that  the  gain  parameters  0^  should  be  zero. 

2.  Computed  Results 

a.  Case  1 

For  this  case  the  assumption  was  made  that  a 
single  component  of  zero  Doppler  and  unit  level  occurred  in 
exact'time  registry  with  the  reference  signal  Rgtt).  Fig.  6 
presents  the  computed  results. 

Because  of  the  time  coincidence  between  the 
signal  and  the  reference  function  R^tt),  a true  null  condition 
requires  that  aQ  be  unity  and  all  other  gain  parameters  be  zero. 
The  manner  in  which  the  condition  is  approached  is  illustrated 
in  Fig.  6 . The  top  graph  in  Fig.  6 shows  the  response  of 
aQ,  0O  as  a function  of  time.  It  is  perhaps  well  to  stress 
that  the  inflections  in  these  curves  are  not  the  result  of  a 
trembling  hand  in  the  graphing  proaess  but  rather  the  result  of 
interaction  among  the  loops. 

The  nominal  behavior  of  the  system  appears  to 

be  as  follows: 

(1)  Because  RQ  is  in  time  coincidence  with 
the  signal,  the  product  of  the  two  signals  includes  a d.c.  term, 
which  on  being  integrated  causes  cxq  to  rise  in  the  positive 
direction,  and  to  introduce  a nulling  signal  into  the  error 
channel.  However,  before  the  null  can  be  completed  the  error 
channel  still  retains  a substantial  level  of  the  received  signal 
and  the  data  suggest  that  during  this  time  the  three  remaining 
channels  act  to  participate  in  the  nulling  action. 
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In  the  description  of  the  Quadrature  Loop 
System  (Section  II-B)  it  was  pointed  out  that  in  an  idealized 

system  a loop  that  is  out  of  time  registry  with  an  FM  signal  by 
a time  € could  still  act  to  null  the  signal  by  describing 
sinusoidal  variations  of  frequency  For  the  case  at  hand 


2 


2tt( 


500 

.128 


) sec*"2 


(44) 


and  c is  a multiple  of  10  msec.  Thus,  the  required  frequency 
of  variation  for  e^.Ol  sec  is 


^ «•*  78.2  Hz 


(45) 


corresponding  to  a period  of  12.8  msec.  In  a nominal  way,  the 
variations  of  and  3^  suggest  the  incorporation  of  frequency 
components  of  this  period  in  their  makeup.  Moreover,  the 
frequency  of  variation  appears  to  increase  appropriately  for  a2, 

P2  and  0,3#  83  to  reflect  their  respective  degrees  of  misalign- 
ment with  the  received  signal.  An  examination  of  Fig.  6 makes 
clear  that  the  gain  parameters  for  other  than  the  zeroth  loop 
cannot  be  represented  exclusively  as  sinusoids  or  even  as  sinusoids 
multiplied  by  a slowly  decreasing  function  such  as  an  exponential. 
The  parameter  3^  appears  to  be  offset  in  the  positive  direction, 
while  a2,  $2»  anc*  £3  are  offset  in  the  negative  direction.  A 
simple  explanation  of  the  noted  offsets  is  not  possible  at  the 
present  time.  However,  the  offsets  appear  to  consist  of  functions 
of  relatively  slow  variation.  One  is  tempted  to  speculate  that 


the  slowly  varying  offset  quantities  could  admit  sufficient  jj 
levels  of  the  various  reference  functions  into  the  error  channel  I 
to  account  for  the  inflections  shown  by  aQ  and  $ 
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For  this  case  the  received  signal  was  assumed  to 
consist  of  three  components,  with  amplitude  distributions  of  8, 

4,  2 and  with  these  components  in  exact  time  registry  with  the 
reference  functions  RQ(t),  R^t),  and  R2(t)  respectively. 

Doppler  was  assumed  to  be  zero  and  the  integrator  time  constant, 
T,  was  assumed  to  be  4 msec. 

In  this  case  again  the  system  converged  to  a valid 
solution,  and  the  values  of  the  gain  parameters  versus  time  are' 
indicated  in  Fig.  7 . Interaction  among  the  loops  is  again  sug- 
gested  by  the  data  and  this  is  again  attributed  to  a tendency  for 
each  loop  to  respond  to  all  of  the  components  of  the  received 
signal.  In  general,  the  parameters  achieve  stable  and  accurate 
values  only  after  the  error  signal  has  been  tightly  nulled. 

c.  Case  3 

The  received  signal  for  this  case  was  assumed  to 
consist  of  a single  component  of  unit  level  whose  onset  time 
occurred  5 msec  after  the  onset  of  RQ(t)  and  5 msec  before  the 
onset  of  R^(t).  Doppler  was  assumed  to  be  zero  and  the  in- 
tegrator time  constant  T was  taken  equal  to  4 msec.  The  response 
of  the  system  to  this  type  of  signal  is  represented  in  Fig.  8. 

In  the  case  under  consideration,  the  system  would 
not  null  the  error  signal.  The  failure  of  the  system  to  achieve 
a null  is  attributed  to  phase  shifts  introduced  into  the  gain 
parameter  signals.  The  phase  shifts  in  turn  appear  to  be  caused 
by  the  integrators  in  conjunction  with  the  feedback  paths  in 
each  loop. 

The  signal  supplied  to  the  error  channel  by  the 
initial  quadrature  loop  can  be  expressed  in  the  form 
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CQ  « a0  cos[wn(t+T5)  + ^(t+T,)2] 


+ 0o  a.tn[uj0(t+T2)  + §(t+T2)2]  . 


The  Bingle  component  signal  is  specified  by 


SQ  « cos[(«0(t+T2-e)  + ^(t+T2~c)2]  , 


where  em5  msec. 


Finally,  the  signal  supplied  by  the  second  system 


loop  is 


Cx  - a x cos[a)0(t+T2-2€)  + ^(t+T2“2e)2] 


+ 0L  sinCuiQ(t+T2-2e)  + ^(t+T2-2e)2]  . 


In  Eqs.  (46)  through  (48),  T2  is  an  onset  delay  that  is  introduced 
to  permit  all  reference  signals  to  assume  their  active  states,  and 
for  the  present  case  ^"30  msec.  The  arguments  on  the  right  hand 
sides  of  Eqs.  (46)  and  (48)  can  be  brought  into  coincidence  with 
that  of  Eq.  (47)  if  the  gain  parameters  are  assumed  to  take  the 
form 
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ao 

Po 

ai 

Pi 


A0(t) 

A0(t) 

Ax(t) 

A^t) 


cos[^e  (t+T2) 
sin[^ic  (t+T2) 
cob[^c (t+T2) 
sin[|ae  (t+T2) 


+ U)q€  - ■Uf-] 

+ <V  " 

3 2 

+ w0e  - -jjie  ] 

+ u)0c  - |ue2]  t 


(49) 


In  Eqs.  (49),  AQ(t) , A^(t)  are  assumed  to  be  constants  or  slowly 
varying  quantities.  Thus,  the  gain  parameters  associated  with 
the  first  two  loops  vary  as  phase  shifted  Binusoids  at  a radian 
frequency  equal  to  pc . 

In  the  case  under  consideration 

w 

^ ■ 2tt  • Tjp  , (50) 

3 

where  is  500  Hz  and  T^  is  128  msec.  Thus,  for  e“5  msec,  the 
frequency  associated  with  the  gain  parameters  is 


JA£ 

2tt 


128 


x 10J  Hz  « 39.0625  Hz. 


(Si) 


The  period  associated  with  thiB  frequency  is  25.6  msec  and  a 
careful  examination  of  the  data  over  time  intervals  removed  from 
the  origin  indicate  that  the  gain  parameters  for  the  first  two 
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loops  vary  at  a frequency  of  39.0625  Hz.  Moreover,  after  the 
parameters  assume  a substantially  sinusoidal  character,  aQ  and 
00  assume  a quadrature  relationship  as  do  and  However, 

the  phase  angles  of  the  waveforms  do  not  agree  with  the  phase 
angles  prescribed  by  EqB.  (49) . For  example,  the  computer 
generated  values  for  aQ  have  zeroes  that  occur  approximately  7.6 
msec  after  the  zeroes  of  the  first  of  Eqs.  (49).  However,  the 
computer  generated  values  for  occur  about  5.2  msec  after  the 
zeroes  of  the  third  of  Eqs.  (49) . The  sum  of  the  two  delays  as 
cited  above  is  one  half  the  period  of  the  sinusoidal  variation 
described  by  the  gain  parameters.  However,  the  significance  that 
should  attach  to  this  is  presently  unknown  and  indeed  the  exact 
causes  of  the  observed  phase  misalignments  are  not  understood. 

The  net  result  of  the  phase  shifts  limitB  the  nulling  action  on 
the  error  signal  such  that  the  rectified  and  averaged  error  sig- 
nal is  about  one  half  the  value  that  would  be  realized  with  the 
received  signal  alone. 
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d.  Case  4 

The  signal  used  in  this  case  consisted  of  three, 
unit  level  components  in  exact  time  registry  with  the  first  three 
reference  signals,  respectively.  The  component  in  registry  with 
Rq (t)  was  assigned  a Doppler  corresponding  to  - 2 Kts  at  25  kHz. 
The  component  coincident  with  R^(t)  was  assigned  zero  Doppler, 
and  the  component  coincident  with  R2(t)  was  assigned  a Doppler 
corresponding  to  +2  Kts  at  25  kHz.  Doppler  was  computed  as  a 
frequency  shift,  A£,  from  the  equation 


Af  *»  2fQ 


Av 

c 


(52) 


where  f^  was  assumed  to  be  the  center  of  an  operating  band  at 
25  kHz,  av  is  a radial  component  of  velocity,  and  c is  the  pro- 
pagation speed  of  sound  taken  to  be  4950  ft/sec.  Evaluation  of 
Eq,  (52)  will  show  that  with  Avss2  Kts,  Af  becomes  34  Hz.  Because 
the  assumed  center  of  the  operating  band  would  be  heterodyned  to 
the  center  of  the  Echo's  Cancellers  operating  band  at  1 kHz,  +2  Kts 
of  Doppler  would  displace  a spectral  line  normally  at  1 kHz  to 
1.034  kHz,  Within  the  signal  generator  subroutine  of  the  computer 
program,  Doppler  is  assigned  by  scaling  the  time  parameter  within 
the  program.  Thus,  for  +2  Kts  of  Doppler,  time  is  scaled  by 
1.034  and  for  -2  Kts  of  Doppler,  time  is  scaled  by  the  factor 
0.966.  The  response  of  the  system  to  the  three  component  signal 
is  represented  in  Pig.  9. 

The  reference  signals  associated  with  a^,  0Q,  and  a2. 

0 2 are  in  time  coincidence,  at  least  initially,  with  the  two 
Dopplered  components  of  the  received  signal.  During  earlier  portions 
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Fig.  9.  Gain  Parameters  Versus  Time  for  a Four  Channel,  Quadrature 
Loop  System.  Case  4. 


of  the  run,  these  parameters  describe  relatively  wide  ranging 
variations  in  response  to  the  error  signal,  but  as  time  progresses 
the  nominal  amplitude  excursions  of  the  parameters  becomes  pro- 
gressively smaller,  and  toward  the  end  of  the  run,  the  small 
excursions  of  these  parameters  indicates  that  correspondingly 
small  levels  of  their  associated  reference  functions  are  being 
admitted  to  the  error  signal  channel. 

The  reference  function  R^(t)  is  in  exact  time  registry 
with  the  second,  zero-Doppler  component  of  the  signal.  The  gain 
parameter,  can  be  seen  to  rise  from  zero  and  increase  toward 
positive  values.  Ideally,  this  parameter  would  approach  unity  to 
null  the  second  signal  component.  However,  as  can  be  seen  in 
Pig.  9 , overshoots  the  value  unity  and  thereafter  describes 
excursions  about  unity  in  response  to  the  composite  error  signal. 
Thus,  can  be  viewed  as  having  developed  a bias  toward  the  value 
that  causes  a null  of  the  second  signal  component,  and  super- 
imposed on  this  bias  are  additive  variations  representing  a response 
to  the  non-zero  components  retained  in  the  error  signal. 
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The  cases  considered  indicate  that  each  of  the  gain 
parameters  shows  some  degree  of  response  to  the  non-zero  com- 
ponents of  the  error  signal.  The  study  has  been  restricted  to 
cases  involving  FM  signals.  For  this  type  of  signal  and  for 
cases  of  zero  Doppler  one  can  also  say  that  each  loop  shows  a 
preferential  response  to  those  signal  components  in  closest  time 
registry  with  its  associated  reference  functions.  The  basis  of 
the  preferential  response  appears  to  reside  in  the  fact  that  the 
gain  parameters  must  describe  a sinusoidal  variation  to  null  any 
component  with  mismatched  time  registry,  that  the  frequency  of  the 
variation  increases  directly  with  the  degree  of  error  in  time 
coincidence,  and  that  each  loop  responds  preferentially  to  lower 
frequencies.  There  exists  a basis  for  considering  each  loop  to 
have  the  frequency  response  of  a single  pole,  low  pasB  filter  with 
a time  constant  of  2T,  where  T is  the  time  constant  of  the  system's 
integrators.  For  a four  msec  time  constant,  the  corresponding 
cut-off  frequency  of  the  low  pass  filter  would  be  about  19.9  Hz. 

In  the  case  Involving  a single  component  mismatched 
by  5 msec  from  "nearest  neighbor"  reference  signals,  the  required 
frequency  of  gain  parameter  variation  was  in  excess  of  39  Hz. 

This  frequency  is  also  above  the  cut-off  frequency  of  19.9  Hz,  and 
difficulties  involving  phase  angle  were  encountered  in  the  case 
in  question.  However,  the  behavior  of  the  system  is  not  yet 
understood  in  a quantitative  sense. 

In  cases  involving  Doppler,  the  system  behavior  becomes 
more  complex,  because  in  this  case  the  frequency  requirements 
imposed  by  a temporal  mismatch  between  a signal  component  and  a 
given  reference  function  can  be  either  increased  or  reduced  by 
the  Doppler  associated  with  the  signal  component. 
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III.  PROPOSED  WORK  FOR  THE  NEXT  QUARTER 


Work  during  the  next  quarter  will  include  a consideration 
of  additional  cases  involving  a reduction  in  the  system's  delay 
line  parameter,  t.  The  balance  of  the  effort  during  the  next 
quarter  will  be  devoted  to  the  preparation  of  a final  report. 


